input(s) to a system and the resultant 
output(s) in real time or a posteriori, or 
from software-generated data sets, were 
presented to the system, which gener- 
ated outputs. Once a system is learned, 
the coefficients and constants can be 


frozen and the algorithm embedded in 
an application. 

This work was done by Michael J. Krasowski 
and Norman F. Prokop of Glenn Research Cen- 
ter. Further information is contained in a 
TSP (see page 1). 


Inquiries concerning rights for the commer- 
cial use of this invention should be addressed to 
NASA Glenn Research Center, Innovative 
Partnerships Office, Attn: Steven Fedor, Mail 
Stop 4-8, 21 000 Brookpark Road, Cleveland, 
Ohio 44135. Refer to LEW-18887-1. 


• Kalman Filter Input Processor for Boresight Calibration 

The new software brings this technology to the industrial level. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Ka-band ranging provides the phase 
center (PC) to phase center range, 
which needs to be converted to the cen- 
ter of mass (CM) to center of mass 
range. Nominally, both PC and CM lie 
on the line connecting the spacecraft 
GRAIL A and GRAIL B. In this case, the 
conversion should be done simply by 
adding the CM-to-PC distance L to the 
measured range for both spacecraft. 
However, due to various technical rea- 
sons, such as displacement of the true 
CM from its nominal position in the 
SRF, or spacecraft attitude fluctuations, 
the PC and CM define a unit vector that 
may be different from the nominal line 
of sight. The objectives of the software 
are to determine the actual line of sight 
direction for each spacecraft and correct 
the previously recorded range data, and 
to provide instructions for how to ma- 


neuver each spacecraft to make neces- 
sary attitude corrections. 

While elements of this approach have 
been used for the boresight calibration 
in the GRACE project, the new software 
brings this technology to the industrial 
level. It is now fully documented and can 
be used by people other than its devel- 
opers. This innovation provides graphic 
outputs and log files that are critical for 
quick analysis and troubleshooting. In 
addition to the line of sight direction, 
the software allows one to evaluate the 
CM-PC base length, which is important 
when the PM location is subject to varia- 
tions (e.g., due to fuel depletion). 

This software is implemented in 
Python and offers excellent cross-plat- 
form porting possibilities. It is very versa- 
tile, and may be applied under various 
circumstances and for other related pur- 


poses. This innovation is capable of com- 
bining the input data from several cali- 
bration maneuvers, evaluating individ- 
ual range biases, and compressing the 
time stamps. It uses Lagrange interpola- 
tion for the orbit data, and a unique 
quaternion-interpolating algorithm for 
interpolating the attitude data. As a re- 
sult, data files with different data rates 
and independent time stamps can be 
handled together. 

This work was done by Dmitry V. Strekalov, 
Gerhard L. Kruizinga, Meegyeong Paik, Dah- 
Ning Yuan, and Sami W. Asmar of Caltech for 
NASA’s Jet Propulsion Laboratory. For more in- 
formation, please contact Brian Morrison at 
Brian.A.Morrison@jpl.nasa.gov. 

This software is available for commercial li- 
censing. Please contact Dan Broderick at 
Daniel.F.Broderick@jpl.nasa.gov. Refer to 
NPO-48479. 


Q Organizing Compression of Hyperspectral Imagery to Allow 
Efficient Parallel Decompression 

Higher compression factors can be attained. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


A family of schemes has been devised 
for organizing the output of an algo- 
rithm for predictive data compression of 
hyperspectral imagery so as to allow effi- 
cient parallelization in both the com- 
pressor and decompressor. In these 
schemes, the compressor performs a 
number of iterations, during each of 
which a portion of the data is com- 
pressed via parallel threads operating on 
independent portions of the data. The 
general idea is that for each iteration it is 
predetermined how much compressed 
data will be produced from each thread. 

A simple version of this technique is 
applicable when the image is divided 
into “pieces” that are compressed inde- 


pendently. As an example, for a com- 
pressor that does not make use of inter- 
band correlation, a piece could be de- 
fined to be an individual spectral band, 
or a fixed number of bands. In the tech- 
nique, the compressed output for a 
piece is comprised of multiple “chunks.” 
The concatenated chunks for a given 
piece form the compressed output for 
the piece. Most of the compressed image 
is produced in multiple iterations, where 
during a given iteration, one chunk is 
produced for each piece. Prior to the 
start of an iteration, chunk sizes are cal- 
culated for each piece. The chunks can 
be produced or decompressed in paral- 
lel. It is noted that it is not specified how 


much of the image data will go into a 
chunk, and in fact a chunk may contain 
incomplete portions of encoded sam- 
ples (at the chunk’s start or end). The 
compressor iterates the process of decid- 
ing on chunk sizes and producing 
chunks for each piece of the requested 
size, until compression of each piece is 
almost finished. At that point, the re- 
mainder of the pieces is compressed se- 
rially without a target chunk size. 

Typically, the chunk size calculation 
should seek to balance the progress 
through each piece, i.e., to leave equal 
numbers of samples remaining in each 
piece; a suggested procedure has this 
aim. A key requirement on the chunk 
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